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Inclusive charmonium production at PANDA experiment 
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The production of the charmonium states in pp experiments is considered at the 
energy rates near threshold in NLO. Such a consideration allows one to obtain non 
zero distributions over the transverse momentum of the final charmonium and gives 

O ' 

■ a natural explanation to the existence of Xci-meson in final state, that is observed 

experimentally and cannot be produced in leading order processes. The crucial role 
of scale parameter choice in a s (Q 2 ) and partonic distributions f(x,Q 2 ) shown, and 
the correct choice offered. 

q_| PACS numbers: 13.75.-n, scattering (energy j 10 GeV) 13.60.Le, 14.40.Pq 

^ : I. INTRODUCTION 

>: 
ov 

Precision experimental studies of the charmonium production in proton-antiproton colli- 
OV n 

sions at low energies proposed in new experiment, labeled PANDA Such studies allows to 
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obtain a deeper understanding of charmonium physics. PANDA experiment deals with pro- 
ton target and antiproton beam energies up to 15GeV. This energies lies near charmonium 
production threshold. In present article we give a theoretical predictions of the inclusive 
charmonium production in pp collisions by accounting next to leading order diagrams in 
partonic cross sections. Such an approach was considered by a number of authors [2, l8r4T0 ] 
and applied mainly for proton-proton collisions. 

The important part of the Program is the detailed analysis of all possible mechanisms of 
charmonia production. Such an analysis is especially important since at low energies there 
is significant difference between charmonium production in pp and pp. For the pp at low 
energies the contributions of gluon-gluon, quark-gluon and quark-antiquark subprocesses 
are comparable, while for pp the quark-antiquark annihilation subprocess dominates. For 
example, if the energy of the proton beam is equal to 40 GeV, the ratio of i/j production 
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cross sections in pp and pp collisions equals &{pp) / o{pp) ~ 6. 

Another problem is that the direct production of ?/>-Hieson is suppressed in comparison 

with the production of the intermediate P-wave states Xco 1 2 with the subsequent decay 

il 

Xc J/ipl- This fact is well confirmed in the experiments [3|. The experimentally observed 
cross sections of Xc2 and Xci are comparable (x c cr meson can hardly be observed due to its 
small radiative width), while the well known Landau- Yang theorem forbids the formation 
of the axial meson from two massless gluons. One more difficulty is that the partonic 
distributions are integrated over the transverse momentum. As a result, such method does 
not allow one to obtain the distributions of Xco- and Xc2-mesons over pt- 

Initially this problems were solved by the introduction of color-octet (CO) components 
of the quarkonia, that arise naturally in the non-relativistic quantum chromodynamics 
(NRQCD), where the expansion over the relative velocity of quarks in meson is performed. 
In this model it is assumed, that final meson is formed from heavy quark pair in color- 
octet state, that subsequently transforms into a physically observed colorless meson. In 
the framework of NRQCD the probabilities of these transitions are described by the matrix 
elements of four-fermion operators, that are determined from the experimental distributions 
over the transverse momentum of final charmonium. We would like to stress, however, that 
this explanation will not work for charmonium production at lower energies. The reason 
is that the distributions caused by octet components decreases slowly with the rise of the 
transverse momentum, but the probability to find such a component in the meson is small, 
compared with the singlet case. As a result, in the large transverse momentum region the 
contribution of octet components can be significant, but for small energies and transverse 
momenta it is suppressed. 

Recently another way to solve this problem was proposed, where the so called non- 
integrated over the transverse momentum distribution functions G(x, h?) are used (/cy- 
factorization) 13Ml5l| . In this case both mentioned above problems are solved simultaneously 



. The transverse momentum of the produced in gluon fusion Xco,2 mesons is explained by the 
transverse momenta of the initial partons. Axial charmonium meson can also be produced 
in gluon fusion, since in the framework of kx factorizations gluons have non-zero virtuality 
of the order k T . There is a number of works, that explain the experimental distributions 
on TEVATRON with the help of these functions (see, for example Il7||). According 



to these works, there is no need to introduce CO components to reproduce the experi- 
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mental data on P-wave charmonium p? distributions. Thus, in /^-factorization approach 
color-singlet components give the main contribution. 

Unfortunately, the method, used in the modeling of the unintegrated distribution func- 
tions G(x, kj<), is based on the summation of large log(l/a;), so it is not applicable for low 
energies, where the gluon momentum fractions are in the range 0.1 < x g < 0.5. For this 
reason we are forced to use the following approximation in our calculations. We start from 
the collinear gluon distributions with well known collinear distribution functions. Further 
we consider the charmonia production at next to leading (NLO) order in the strong coupling 
constant a s . Such a trick enables us to obtain the distributions over p? for all charmonium 
states. For Xco and Xc2 production we observe a collinear singularity at px = 0. To avoid 
this singularity we introduce a regularization procedure. For directly produced ip and Xci 
such a singularity is absent and we use the whole integration region for p?. 

Next section is devoted to the consideration of different modes of charmonia production 
and collinear singularities regularization. In the third section we determine the cross sections 
of the hadronic processes and pay attention to the correct scale parameters choice in a s (Q 2 ) 
and parton distributions f(x,Q 2 ). There we are present numerical results. 



II. PARTONIC SUBPROCESSES 



Feynman diagrams corresponding to the charmonium production are presented in FigJTJ 
At the leading (FigHK) order only processes gg —> Xco,2 are available. Process gg — >■ if) is 
forbidden by charge parity. Process gg — > %& is forbidden due to Landau- Yang theorem, 
that forbids the formation of the axial meson from two massless gluons. Other LO cross 
sections we give in Appendix lAl 

Next to leading order diagrams are presented in FigJTb,c. First two diagrams in FigJIL, 
includes 3-gluon vertex. To avoid delicate problems, bounded with ghosts contributions, we 
recalculated differential cross sections for process gg —> Qg (Q = if), Xeo,i,2) (and also for 
qg -> Xo,i,2?) in axial gauge: 

Cg f = -^(n,A a ») 2 , with n/ = -1 (1) 

This gauge does not require additional ghosts Lagrangian, but gluon propagator and po- 
larization sum becomes more complicated. Using £ = (Landau choice), it can be found 
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Figure 1. Feynman diagrams for the charmonium production, a) direct Xc0,2 production, b) 
through qg — > Qq and qq — > Qg subprocesses, where Q = XcJ- c) through gg — > Qg subprocess, 
where Q = XcJ for the first two diagrams and Q = ip, XcJ for the last diagram. 



that 



k^k u k^riy + k u n^ 



(k,n) 



(k,n) 
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Auxiliary vector n„ will disappear in physical observables. We found, that our results are 



in excellent agreement with 9|, 10| for gg channel and |2j for qg channel. Exact formulas 
for the differential cross sections are rather tedious and can be found in the cited papers. 

Massless t-channel gluon in propagators in Fig{T] leads to the collinear singularities in 
small t and u = M 2 — s — t regions, so gg — > Xo,2<7 and qg — > Xo,2l are divergent: 
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dt 
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In terms of meson transverse momentum p T = ytu/s, these singularities correspond to 
Pt — >• singularity. To calculate the total cross section: 



o is) 



(s-M 2 )/2 v / I 



dt 



dp T 



(4) 



M 2 -s 

some regularization should be performed. Popular decision is to cut off the small pt region, 
i.e. restrict integration region in last formula by setting p? > A, where cutoff parameter A 



can be taken from experimental setup or from some physical reasons. For example in [2[, A 
was taken equal to 1/ R c5 , where R C c is the geometrical size of charmonium. Such approach 
have a big lack, because the total cross section is high sensitive to A variation. To avoid 
these difficulties, we will use another regularization procedure. 

Taking indefinite integral of differential cross section, the following well known relation 
can be found: 

f d&{gg^Qg) a s I ' M 2 \ u 

at - 2 = —a (gg ->• Q)P g ^ gg — In T + finite (5) 

at Zn \ s ) t 

Similarly, for qg — > Qq reaction we have 

di — = ^M S9 -» fl)P„ f-) In f + finite. (6) 

dt 2n \ s / t 

In these expressions the second parts are finite at t — > and t — > M 2 — s, while &o{gg — > Q) 
are given in (1A2[) . f1A2h . P g ^ gg (x) and P q ^ qg (x) are well known QCD splitting functions: 

(x 1 — x \ 
1 h x(l - x) 
1 — X X J 

p (x) - ll+il^l! 

The full hadronic cross section can be obtained by integrating partonic cross sections with 
partonic distribution functions: 

a(s) = J dx 1 dx 2 f(x 1 )f(x 2 )d-(s) (7) 

Singular parts of ()5]) and ([6]) are included in partonic distributions functions f(x) and gen- 
erate well known scaling violations, described by Altarelli-Parisi equations. So, inclusion of 
singular parts in partonic cross sections leads to double-counting: one time in a and one in 
f(x). Thus, we will use regularization: 

do{gg ->• Qg) z a s ( M 2 \ w N " 



-«»)=(/ mam^m di _ ^ {gg _> Q)P ^ g ^ j lB 

and similarly: 



(8) 

t=M 2 -s 



t=0 

(9) 

t=M 2 -s 



In contrast to Xco,2, differential cross sections for if) and Xci have no collinear singularities 
and finite at t — >• and u — >■ 0. In case of Xi, this is explained by the Landau- Yang theorem, 
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that forbids the production from two massless gluons. As a result, the squared matrix 
element of this reaction is proportional to the virtuality of the intermediate t-channel gluon, 
so this factor compensates the divergency, caused by the propagator. For ip, we have similar 
reasoning based on charge parity, since the first two diagrams of FigJTb are absent in this 
case. 

Exact formulas for regularized partonic cross sections are given in Appendix [A] 



III. HADRONIC CROSS SECTIONS 

Let us now consider full hadronic process 

A(P 1 )B(P 2 )^Q(P)+X, (10) 

where A and B are the initial hadrons, Q = ip, Xcj , and in the parentheses corresponding 
particle momenta introduced. The cross section of this reaction is expressed through the 
cross sections of considered above partonic reactions : 

a ( S ) = Yl J dx 1 dx 2 fa/A(x 1 )f a / B (x2)o-ab(s), (11) 
a,b 

where summation is performed over partons a and b, x± }2 are the momentum fractions held 
by these partons, and f a /A{xi), fb/B(x 2 ) are the distribution functions of the partons in the 
initial hadrons. In the common variables 

x = X\ — x 2 (12) 

s = (x 1 P l + x 2 P 2 ) 2 = xix 2 s, (13) 



the full hadronic cross-section becomes 

s x(s) 

1 1 r 

X = Xi+X 2 , (14) 



a ' b M 2 -x(s) 



Xl,2=Xl } 2(x,s) 



where 



x(s) = l-- (15) 
s 

In our numerical estimates we used the distribution functions and a s numeric values pre- 



sented in the work 



jll |. Other numerical parameters are equal to: 

= 3.097 GeV, M Xc0 = 3.415 GeV, (16) 

M Xcl = 3.511 GeV, M Xc2 = 3.556 GeV (17) 

#1(0) = 0.81 GeV 3 , #p(0) = 0.075 GeV 5 (18) 
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Figure 2. a) Q-dependence of u-quark distribution at x = 0.9 b) Q-dependence of a s 

A. Scale dependence 

There are two physical quantities, that depend on some scale choice: partons distribu- 
tions f a /A{x,Q 2 ) and strong coupling constant a s (Q 2 ). The parton distribution function 
fa/A(%,Q 2 ), gives the probability of finding a parton a of longitudinal fraction x in phys- 
ical (anti)proton, taking into account collinear gluon (or massless quark) with transverse 
momenta px < Q- It is clear, that exact value of Q depends on parameters of partonic 
subprocess, i.e. Q 2 = Q 2 (s). It is convenient to set Q 2 to a fixed value Q 2 - characteris- 
tic momentum transfer of partonic subprocess. Such choice can be argued by mean value 
theorem, which states: 

b b 

J f(x)g(x)dx = J g(x)dx 

a a 

On the other hand, from the structure of Altarelli-Parisi equations, it is clear, that at least 
for high energies (Q > 1 GeV) error in choosing of leads to negligibly small variation of 
final results; for example, the parton distributions change by ~ 1% as Q* is changed by a 
factor of 10. Another situation we have at low energies (Q < 1 GeV), when perturbation 
theory works bad, and error in choosing of leads to a dramatic variation of partonic 
functions. At the energy rates ~ 1 GeV, error in Q* in two times can leads to f(x, Ql) error 
about 20%. For example, u-quark distribution dependence on is shown on Figj2]a). 

Second physical quantity, depending on some scale choice is a s {Q 2 ). Here, the meaning 
of Q is different from the meaning in distribution functions. The strong coupling a s scale 
dependence occurs when the full propagators and vertexes are inserted in tree level diagrams. 
Similar reasoning, based on mean value theorem allows to set a s {Q 2 ) to a fixed value with 



8 



some scale Q*. The perfect justification of this procedure is given in [12| . Only in simplest 
situations the exact value of can be found: for example for 2—7-2 reaction through 
s-channel particle, it can be found, using Callan-Symanzik renorm-group equation, that 
Q 2 = s. In other cases the exact value of scale can be found only by sequential analysis 
of perturbation series expansion. It is clear, that Q* depends on the process. At high 
energies Q > 1 GeV a s (Q 2 ) becomes almost constant. At the energies rates near 1 GeV, 
a s dependence of Q is shown on Figj2]b). 

So, in general, we have three possible ways of setting scale parameters in full cross section: 

fixed scheme a(s, Ql) = £ a ,6 / ds a ab (s, a s (Ql)) J ff a /A(x 1 ,Ql)f b/B (x 2 ,Ql) 

M 2 -x(s) 
s x(s) 

float scheme cr(s,Q 2 ) = J2 a ,b I ds a ab (s, a s (s)) J ! ff a /A(xi,Ql)f b / B (x 2 ,Q 2 J 

M 2 -x(s) 

float2 scheme a(s,Ql) = E a> 6 / ds a ab (s, a s (s)) / ^f a /A(xi,s)f b/B (x 2 ,s) 

M 2 -x(s) 

The fixed scheme is the most common way, used in calculations. The float scheme, takes 
in account the fact, that a s in partonic subprocess, depends on s. Of course, in general, 
this dependence is complicated and has the form a s (f(s)), but from general considerations 
it is clear, that for small interval, near process threshold, f(s) ~ s. The float2 scheme, 
takes in account both a s (Q 2 ) and f(x, Q 2 ) scaling. It is also clear that maximum transverse 
momentum Q in f(x, Q 2 ) depends on the s and for small energy intervals we set Q 2 ~ s. 
Actually, the cross section in float2 scheme does not depends on Q*. 

Figj3] illustrates the dependence of the cross section on scheme choice for Xco production 
in uu channel at pp collisions. It is seen, that all three curves are crossed in one point at 
Q2 = M 2^ = ii. 66 GeV 2 . For other mesons and other channels the picture is equivalent - 
three curves are crossed at the corresponding squared meson mass. The only exception - 
quark-gluon channel, where cross point greater at 10%, then corresponding squared meson 
mass (FigHJ). We make similar calculations for other energy regions and found, that this 
results remain valid. As was excepted, at high energies (s 3> lGeV 2 ) the difference between 
schemes is negligible. So, this results improves that all schemes are equivalent with the 
correct Q* choice: 

Q* = M, (19) 

where M - corresponding meson mass. 

In all further calculations we shall use the fixed scheme with = M. 
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Figure 3. Scale dependence in different schemes of the uu channel cross section in pp — > XcoX 
collision at the energy ^/s = 4.34GeV. 
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Figure 4. Scale dependence in different schemes of the ug channel cross section in pp — > x&X 
collision at the energy yfs = 4.34GeV. 



B. Total cross sections 

Figj5] illustrates the dependence of if) production through different processes on total 
energy. The bold line shows the summed over all processes cross section: 



° = VggW + Br (XcO -»■ V'tMXco) + Br(x c i -»- H)v(Xci) + Br(x c2 -»■ V^Mx^), (20) 
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o ( s ) , pb 




Figure 5. Different contributions to the total tp production in proton- antiproton reaction for 
different cm. energies. 1 — total ip production, 2 — ip production through the radiative Xc2 
decay, 3 — through the radiative Xd decay, 4 — through the radiative Xco decay, 5 — direct tft. 

where the branching values are equal to: 

Br( Xc0 -> H) = O- 016 . Mxa -> = 0-344, Br( Xc2 -»- H) = 0.195 (21) 

From this picture it is seen, that direct ip production and production from radiative Xco 
decay is highly suppressed. The contribution of Xco radiative decay is negligible, due to 
very small branching value. As was noted above, the direct ip production only available 
in process gg — > tpg, but the gluon-gluon channel in pp reactions is highly suppressed by 
the quark- ant iquark channel, so the cross section of the direct ip production is significantly 
smaller, then the production of Xcj, where the quark-antiquark channel is available. 

In the Fig IH] and Fig 17] we show the contributions of the different subprocesses to the total 
Xd and Xc2 production cross sections. For both mesons the most significant contribution 
gives the uu subprocess. For the gluon-gluon subprocess our numerical results are equals 
to zero within the error of numerical calculations. The negligibly small effect of the other 
channels can be easily explained by the structure of the parton distributions. The small 
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energy of the hadronic reaction corresponds to the longitudinal fraction x in partonic dis- 
tributions fa/A(x) close to zero. For this region the u-quark distribution function absolutely 
dominates. 



C. Production mode at the s = 32 GeV 2 



It should be stressed, that presented above expressions for charmonia production cross 
sections can be considered only as estimates on upper bounds. The reasons is that in some 
events initial (anti)protons are present also in the final state. Due to baryonic number 
conservation in proton-proton scattering this configuration is realized almost always. In 
proton-antiproton interaction, however, that presence of baryons in the final state is not 
necessary. Numerically this effect can be described in terms of inelasticity coefficient, which 
can be interpreted as the probability of proton-antiproton annihilation into other states. 
According to presented in [18] analysis, this coefficient is equal to K ~ 0.5 and decreases 
slightly with the increase of energy. 

If initial baryons are present also in the final state, the effective interaction energy de- 
creases from y/s to ^/stf{ ~ \fs — 2M p . In the case of high-energy colliders this modification 
does not change significantly the cross sections of the considered processes. For PANDA en- 
vironment, however, situation is completely different. From fig|5]it is clear, that the decrease 
from y/s ~ 5.5 Gev to y/s c s ~ 3.5 GeV leads to dramatic decrease of charmonia production 
cross sections. So, one can expect, that the reactions pp — > pp + J /if) + X give negligible 
contributions to the cross sections of charmonia production at PANDA, and expression ( I14p 
should be multiplied by the inelasticity factor K ~ 0.5. 

At the production mode in PANDA experiment the antiproton beam energy equals to 
15 GeV, that corresponds to the s value equals to 32 GeV 2 . The cross section of if) meson 
production is given in ( 1201) . Our calculations give 

a(pp^ ipX) = 0.21 nb, (22) 

where 



a{pp-+ XdX) = 0.2 nb 
a(pp — > Xc2X) = 0.75 nb 
a(pp XcoX) = 0.35 nb. 
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o ( s ) , pb 




Figure 6. Contribution of the different subprocesses to the total Xci production. 1 — uu subprocess, 
2 — dd subprocess, 3 — the sum over ug and ug subprocesses, 4 — the sum over dg and dg 
subprocesses. 

o ( s ) , pb 




Figure 7. Contribution of the different subprocesses to the total Xc2 production. 1 — uu subprocess, 
2 — direct Xc2, 3 — dd subprocess, 4 — the sum over ug and ug subprocesses, 5 — the sum over 
dg and dg subprocesses. 
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The ratio of Xci and Xc2 production cross sections is equal to 



<?{x 



0.26 



c2 



(23) 



The pt distribution of mesons, can be obtained by rewriting differential cross sections in 
terms of Pt = Jiu/s and integrating with the partonic distributions: 



da 
dp T 



where 



(Pt + a/pt+M 2 ) 2 



da 



ds da(ab — > Qc) 
s dpx 



2sp 7 



^-fa/A(x 1 )f b/B (x 2 ) 



(24) 



-(i-f) 



dp T 

tl,2 



s - M 2 ) 2 - Asp 



/da 




da 




\ di 


t=ti 


di 


i=t 2 / 



-(M 2 -§±yj(§- M 2 ) 2 - Asp 2 T ) 



As was shown before, the major processes, giving contribution to the total ip production, 
are the radiative decays of the Xa,2 mesons, which in turn are formed through the uu 
subprocess. Thus, for the calculation of the pt dependence we will neglect all other channels 
of the ip production. In thus approximation we do not encounter with collinear singularities, 
that appear in other subprocesses of Xci formation. Another problem arises when we consider 
the total ip distribution. The radiative decays Xcj ^7 can gi ye significant contribution to 
the ]>r-distribution of the final if), when the transverse momentum ~ 1 GeV. However, we 
will neglect such contribution. On the FigJH]we show the transverse momentum distributions 
of the if) production: 



dp T 
da(ip) 

dp T 



ds da(uu ->• Xa,2g) f dx 



dpi 



X 



c i ^ i \ da (Xd) , uf _. ,da{xc2) 



dpq 



dpi 



IV. CONCLUSIONS 



The paper is devoted to J /if}- meson production in proton-antiproton interaction at low 
energies. This process can be used to clarify modes of charmonia production in hadronic 
experiments and allows one to measure with higher accuracy proton spectral functions at 
x ~ 0.5. 
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1 da nb 
p T dp r ' GeV 2 




Figure 8. Transverse momentum distributions of ip production with inelesticity coefficient taken 
into account. 1 — total ip, 2 — ip production through Xc2 decay, 3 — ip production through Xd 
decay 

The physics of charmonia production in hadronic reactions is completely different for 
different energies. For high-energy experiments (e.g. Tevatron or LHC) heavy quarkonia 
are produced mainly in the gluon-gluon interaction, since small values of feynman variable 
x are allowed kinematically. The contributions of quark-gluon or quark-antiquark modes 
are negligible. In the threshold region, where only values x ~ 0.5 are allowed, on the 
contrary, main contributions come from quark-gluon (in proton-proton interaction) or quark- 
antiquark (in the proton-antiproton interactions). In the nearest future at the FAIR proton- 
antiproton particle accelerator with 3 < y/s < 5.7 GeV the PANDA detector will perform 
first measurements, so a reliable prediction for charmonium mesons production for this 
experiment is required. 

In our paper we give predictions for total cross sections of J/ip-meson production in 
different modes at NLO. The emission of additional gluon leads to non-zero transverse mo- 
mentum of final charmonium, that is obviously absent in LO partonic reactions gg — > Xco,2- 
It is shown, that main contributions to this process are given by Xd^-mesons production 
due to quark-antiquark annihilation with the subsequent radiative decay Xci,2 — > J 1^1- 
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Special attention is given to regularization of infrared and collinear singularities in the 
case of Xc2 meson production, when the t-channel gluon in gg — > Xc29 partonic reactions 
leads to divergency in pt distribution and infinitive values of the cross section. 
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In this section we give the total cross sections of partonic subprocesses. For subprocesses, 
in which the collinear singularities appears, we use the regularization procedure, described 
in the main text. 



In the leading order only gg — > Xco,2 reactions are possible. The cross sections of these 
reactions are 
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Appendix A: Regularized partonic cross sections 



1. Leading order 



(Al) 



where 



vo(gg -»■ xo) = 12 



7r 2 a 2 R' 2 x (0) 



M 5 s 



<to(gg ->■ X2) = 16 



7r 2 a 2 R' 2 (0) 



M 5 s 



where R' x (0) is the derivative of the radial part of x- meson wave function at the origin, 
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2. gg -> 

As was noted above, the cross sections for processes gg — > ipg and gg — > Xci9 have no 
collinear singularities, so they does not require regularization/ We have 

where -R</,(0) is the radial part of ip wave function at the origin, and 



£(00 -> = * ( 12M 4 s (M + 9M 14 s+ 

26M 12 s 2 + 28M 10 s 3 + 17M 8 s 4 + 7M 6 s 5 - 40M 4 s 6 - 4M 2 s 7 - 4s 8 ) log 

J s 

(M 2 - s)(M 2 + s) (M 18 + 39M 16 s + 145M 14 s 2 + 251M 12 s 3 + 119M 10 s 4 - 
153M 8 s 5 - 17M 6 s 6 - 147M 4 s 7 - 8M 2 s 8 + 10s 9 



The gg production of Xco,2 states, have collinear singularities. Performing regularization 
procedure (|S]), we obtain 

&**(ag -> Xo ^) = 27m 3 i^(0) , i32M 22 s - 7M 20 s 2 - 

vyy Aoy; 3M 7 s 3 (M 2 — s) (M 2 + s) V 

80M 18 s 3 + 210M 16 s 4 - 560M 14 s 5 + 802M 12 s 6 + 696M 10 s 7 - 

1721M 8 s 8 - 244M 6 s 9 + 789M 4 s 10 + 56M 2 s n + 12s(-24M 22 

-41M 20 s + 10M 18 s 2 + 7M 16 s 3 + 42M 14 s 4 - 176M 12 s 5 - 10M 10 s 6 + 

M 2 \ 

40M 8 s 7 + 14M 6 s 8 - 31M 4 s 9 + 9s 11 ) log — - 172s 12 J , 

a Re9 (gg -> XzQ) = f (66M 24 + 201M 22 s- 

vyy Aiy; 3M 7 s 3 (M 2 — s) (M 2 + s) V 

31M 20 s 2 - 728M 18 s 3 + 360M 16 s 4 - 266M 14 s 5 + 256M 12 s 6 + 1032M 10 s 7 - 
752M 8 s 8 - 271M 6 s 9 + 207M 4 s 10 + 32M 2 s n - 12s (l2M 22 + 5M 20 s + 17M 18 s 2 + 
86M 16 s 3 - 204M 14 s 4 + HM 12 s 5 + 31M 10 s 6 + 74M 8 s 7 - 8M 6 s 8 + 22M 4 s 9 - 
Qs u j log — - 106s 12 J . 
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3. qg ->• 

At the qg channel the only Xcj mesons can be produced. The Xi meson cross section 
does not requires regularization and is equals to 

x 167ra 3 i^(0) / . „ „ , s ,\ . . 

HQ9 -+ Xiq) = QM7g r ( 4M - 9M § + 3M § lo § M2 + 5S J • ( A2 ) 

Regularized cross sections for Xco,2 are 

^ 9 (q S Xo<?) = - 16 ;g^ (0) (4M 6 - 18M 4 s + 57M 2 s 2 + 
3s (4M 4 - 9M 2 s + 9s 2 ) log - 43s 3 ) , 

27M 7 s 3 

3s (5M 4 - 12M 2 s + 12s 2 ) log - 53s 3 ) . 



* Re9 (qg X 2 g) = - 16 :"1^ (0) (20M 6 - 36M 4 s + 69M 2 s 2 + 



4. -»• Qs> 

At the channel, all cross sections are finite. This is explained by the fact, that their 
differential cross sections are cross-symmetric (t -H- s) to the qg ones: 

\M(qg^Qq)\ 2 = \M(qq ^ Qg)\ 2 \^ s (A3) 

and the total cross sections are: 

1287ra 3 i? /2 m) , 9X 2 



2567ra 3 J R'V0) , , 9 9 \ 

* ( « - **> = - W (tf -i) ( 6M + 3M s + s ) < A5) 

256vra 3 J R ,2 ,(0) , _ N 
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